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A beamformvr is a processor used in conjunction with an 
array of sensors to provide a versatile form of spatial filtering. 
The sensor array collects spatial samples of propagating wave 
fields, which are processed by the faeamformer. The objective 
is to estimate the signal arriving from a desired direction In 
the presence of noise and interfering signals. A beamformer 
performs spatial filtering to separate signals that have over- 
lapping frequency content but originate from different spatial 
locations. This paper provides an overview of beamforming 
from a signal processing perspective, with an emphasis on re- 
cent research. Data independent, statistically optimum, adap- 
tive, and partially adaptive beamforming are discussed. 



I, INTRODUCTION 

The term beamforming derives from the tact that early 
spatial filters were designed to form pencil beams (see 
poldT plot in Fig. 1 .1) in order to receive a si^^nal radiating 
from a specific location and attenuate signals from other 
k)calions. "Forinlnji beams" seems to indicate radiation of 
energy; however, beamforming is applicable to either 
radiation or reception of energy. \n this paper we discuss 
formation of beams for reception. 

Systems designed to receive spatially propagating sig- 
nals often encounter the presence of interference signals. 
If the desirod signal and interfcrers occupy the same tem- 
poral frequency band, then temporal filtering cannot be 
used to separate signal from interference. However, the 
desired and interfering signals usually originate from dif- 
ferent spatial locations. This spatial separation can be ex- 
ploited to separate signal from interference using a spatial 
filter at the receiver. Implementing a temporal filter requires 
processing of data collectt?d over a temporal aperture. 
Similarly, implementing a spatial filter requires processing 
of data collected over a spatial aperture. 

Several applications that employ spatial filtering of data 
are listed in Table 1.1. Fig. 1.1 illustrates a microwave com- 
munications antenna that employs a continuous spatial 
aperture to accomplish spatial filtering with a single an- 
tenna, fig. 1.2 depicts a low frequency-towed sonar array 
in which the spatial aperture is obtained through a dis- 
crete spatial sampling by an array of sensors. When the 
spatial sampling Is discrete, the processor that performs 
the spatial filtering Is termed a beamformer. Typically a 




Figure 1.1 A continuous spatial aperture provides one 
mechanism for spatial filt:er<ing. Illustrat;ed is o parabolic 
microwave antenna system. The antenna dish provides 
the spatial aperture over which energy is coiiected. This 
energy is reflected to the antenna feed. The dish and feed 
operate as a spatial integrator. The energy from a far field 
source located directly in franc of the antenna arrives at 
the feed temporarily aligned (i.e.. all source-to-feed path 
lengths are equaD and is coherently summed. In general, 
energy from other sources wilt arrivs at the feed via 
variabia length paths, and add incoharantly. A polar plot 
of a typical antenna beampattem (i.e., power gain vs.. in 
this case, azimuch angle) is shown for a selected fre- 
quency and for the elsvation angle at which the antenna is 
pointed. 
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TABLE 1.1 




ARRAY6 AND BEAMFORMERS PROVIDE AN EFFECTIVE AND VERSATILE MEANS OF SPATIAL 


FILTERING. THIS TABLE USTG A NUMBER OF APPLICATIONS OF SPAHAL FILrEflINQ, GIVES 


EXAMPLES OF ARRAYS AND 8EAMF0RMERS. AND PROVIDES A FEW KEY REFERB^CES. 


Appltcation 


Description 


References 


RADAR 


phased-array RADAR; air traffic 


Brookner (19851; Haykin [1985]; 




control; synthetic aperture RADAR 


Munson et al. {19831 


SONAR 


source localization and classification 


Knight et al. I1981J; Owsley 






ri98Sl 


Com mu ni cations 


directional transmission and * 


Mayhah [1976]; Compton 




reception; sector broadcast in 


[1973]; Adams et al. 11980] 




satellite communications ' 




Imaging 


ultrasonic; optical; tomographic 


Macovski.11983]; Pratt 119781; 




Kak [1985] 


Geophysical 


earth cni5t mapping; oil exploration 


lustlce. [19851 


Exploration 






Astrophysrcal 


high resolution imaging of the 


Readhead [19821; Yen [1985 


Exploration 


universe 




Biomedical 


fetal heart monitoring; tissue 


Widrow et aJ. {197S1; Gee et al. 




hyperthermia; hearing aids 


[1984]; Peterson et al. [1987] 



l^camformer linearly combines the spatially sampled time 
scries irom each sensor to obtain a scalar output time 
scries in the same manner that an FIR filter linearly 
combines temporally sampled data. Two principal advan- 
tdfc;e5 of spatial sampling with an array of sensors are dis- 
cuisod hc'low. 

Spatial discrimination capability depends on the size of 
the spatial aperture; as the aperture fnr.reases, discrimi- 
nation improves. The absolute aperture sizn is not impor- 
tant, rather its size in wavelengths is the critical parameter. 
A single physical antenna (continuous spatial aperture) 
capable of providing the requisite discrimination is often 
practical for high frequency signals since the wavelength 
is short. However, when low frequency signals are of in- 
terest, an array of sensors can often synthesize a much 
larger spatial aperture than that practical with a single 
physical antenna. 

A second very significant advantage of using an array of 
sensors, relevant at any wavelength, is the spatial filtering 
versatility offered by discrete sampling. In many applica- 
tion areas it is necessary to change the spatial filtering 
function in real lime to maintain effective suppression of 
interfering signals. This change is easily implemented in a 
discretely sampled system by changing the way in which 
the beamformer linearly combines the sensor data. 
Changing the spatial filtering function of a continuous ap- 
ifrture antenna is impractical. 

The purpose of this paper is to describe bcamforming 
from a signal processing perspective, provide an overview 
of beamformer design, and briefly discuss perform- 
ance and implementation issues with an emphasis on re- 
cent research. The paper begins with a section devoted to 



defining basic terminology, notation, and concepts. Sue- 
ceeding scc*trons cover data-independent, statistically op- 
timum, adaptive, and partially adaptive beamforming. We 
then provide a brief discussion of implementation issues 
and conclude with a summary. 

Throughout the paper we use familiar methods and 
techniques frorn FIR filtering to provide insight into vari- 
ous aspects of spatial filtering with a beamformer. How- 
ever, in some ways bcamforming differs significantly from 
FiR filtering. For example, in beamforming a source of 
energy has several parameters that can be of interest: 
rangc^ azimuth and elevation angles, polarization, and 
temporal frequency content. Different signals are often 
mutually correlated as a result of muttipath propagation. 
The spatial .sampling is often nonuniform and multi- 
dimensional. Uncertainty must often be included in char- 
acterization of Individual sensor response and location, 
motivating development of robust bcamforming tech- 
niques. These differences indicate that bcamforming rep- 
resents a more general problem than FIR filtering and as a 
result, more general design procedures and processing 
structures are common. 

Rather than making a futile attempt al attributing devel- 
opments due to many different researchers in beam- 
forming, wc refer the reader to the following references: 
books — J. W. R. Griffiths et al., ed. [1973), Hudson 119811, 
Monzingo and Miller [1980], Haykin. ed. (19851, Compton 
(1988]; special issues — tE££ Transactions on Antennas and 
Propagation [1976], (1986), Journai of Ocean Engineering 
[19871; tutorial — Gabriel 119761; and bibliography— Marr 
11986J. Papers devoted to beamforming are often found in 
the IEEE Transactions on: Antennas and Propagation, 
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Figure 1 .2 This figure exempfifiBS^^spsitJal-filtQ^ iS; 
discrete spercure from a towed SOI^lAB; array ot^^ 
spaced sensors. The signalis recej^ed ettto sensof^i^-; 
Co far field sources' locaeed acarr^y.brpaCteidei;:G.e;^:fira 
aiong the plane perpendictiiar to the ierray-«idB)v witti^ 
temporarily aligned across the anSay/assuit^ 
sors are Identical end propagatiah 19 isotr&jpikL^i^n^^ 
thg sensor outputs results : in d coherieht^ jg|an s^ferji^^^^ 
that demonstrated In Figure 1 . 1 . Jf 8 aoOP6eifremflrtfit*^ 
DirectiorHDf-ArnvalL. CDCM) is of; iiftt^resls ife])^^ 
electronicaJly implamented at the .sensor DtIjt;put6L!to:'^^^ 
pensac& for that iDCM^s intdre^ns^ li^rq^diatidft 
i^iereby .temporartly aligning: the aquirpe^. :7Tie/d6llaj^ 
sunvning device constrtutb a bBamforTT«r.:;A -r^ 
satila approachis to usis f jtters at the asns^. oii^EJSift:^^: 
sum the fitter; outputs,: A':r^ctangulacbp).bt::of itJie^'ra 
n'(tud&^^.■DaA•i8 shown fbr a phase shi1!trwjd| sumtia^^ 
former. A linear equi^paced array of ;i^^cMlcfiMat<^^ 
sensors is enipidyed. The be«imforimr Is '^Bt^hwTrta '18^. - ' 



Acoustics. Speech, and Signal Processing, Aerospace and 
Electronic Systems, and in the loumaf of the Acoustical 
Society of America. There is a vast body of literature on 
various aspects of bcamforming and we can only refer- 
ence a subset in support of our discussions. We often 
refer io the FIR filtering literature in our discussions of 
beamforming, since their histories are both parallel and 
overlapping. 

II, BASIC TERMINOLOGY AND CONCEPTS 

In this section we introduce terminology and concepts 
employed throughout the paper. We begin with a subsec- 
tion that defines the beamforming operation and discusses 
spatial fihering. The next subsection, entitled "Second 
Order Statistics." develops representations forthecovari- 
ance of the data received at the array and discusses dts> 



tinctions between narrowband and broadband beam- 
forming- The final subsection defines various types of 
beamformers. 

A. Beamforming and Spatiaf Filtering 

Fig. 2.1 depicts two beamformers. The first, which 
samples the propagating wave field in space, is typically 
used for processing narrowband signals. The output at 
time k, y{k}. is given by a linear combination of the data at 
the J sensors at time k: 



yik) - 2 w,*x.(k) 



<2.1) 



where * represents complex conjugate. It Is conventional 
to multiply the data by conjugates of the weights to sim- 
plify notation. We assume throughout that the data and 
weights are complex since In many applications a quadra- 
ture receiver is used at each sensor to generate in phase 
and quadrature (I and Q) data. Each sensor is assumed to 





f igi^^S^^lr iiAriMam^ combination of 

thesefMvmte^^ are each multi- 

pU^t»y^a«!fBpi(^ then suriOTiBd [for nbtational 

cpiivMKnceJt4s ooft^tion^ to apply conjugates 6f the 
dafiwcln^^i^l^^ Is typically used with 

;i)4mMi(Brtf beamformer is 
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have any necessary receiver electronics and an A/D cori- 
vertcf if beamforming is performed digitaHy. 

The second beamformer in Fig. 2.1 samples the propa- 
gating wave field in both space and time and is often used 
when signals of significant frequency extent (broadband) 
are of interest. The output In this case can be expressed as 



y(k) = i wCpX,(k - p) 

1-1 p-O 



(2.2) 



where K - 1 is the number of delays in each of the J sensor 
channels. If the signal at each sensor is viewed as an Input, 
then a beamformer represents a multi-input Single out- 
put system. 

It is convenient to develop notation which permits us to 
treat both beamformers in Rg. 2,1 simultaneously. Note 
that (2.1) and (2.2) can be written as 



y(k) = w^x{k). 



(2.3) 



by appropriately defining a weight vector w and data vec* 
tor x(k). We use lower and upper case boldface to denote 
vector and matrix quantities respectively, and let super- 
script H represent Hcrmitian (complex conjugate) trans- 
pose. Vectors are assumed to be column vectors. Assume 
that w and x(k) are N dimensional; this Implies that K = KJ 
when referring to (2.2) and N ^ J when referring to (2.1). 
Except for Section V on adaptive algorithms, we vyill drop 
the time index and assume that its presence is understood 
throughout the remainder of the paper Thus (2.3) is writ- 
ten as y = vv^x. Many of the techniques described in this 
paper are applicable to continuous time as well as discrete 
time beamforming. 

The frequency response of an FIR filter with tap weights 
w;;, 0 ^ p ^ J and a tap delay of T seconds is given by 



Alternatively 



r(w) = 2 vvje" 



r(tt>) = w"d(cd) 




(2.4a) i 



(2.4b) 



IfjQtjr^^i^ .^i'^^ )iutthy3ached.delay lines providBS a 

i^fe aifalfiy ^ sources. TtMs 

^tlfljtjc^i^^ a siwwlrpropagating in 

^p^eil^ii^giKt^ra^piBQc^^ at uOA 0. With J serv 

l^^w^i^s^^^ the 

^^^nS|loi^if«<ji^ at :A/ = J/C noh- 

fiii^il^^d^^'pgpp^i rt;^ from the 

,|;|M^^|toBf:it^!^^ of the 

iJ^iWflajrijIi^i^^ apefture of the obsar- 

lH!«6ji6|fepi^^ suggests, tem- 

Li:iji-i=i.-^™iji^^jj^;,^j^^^ 0 pigno wave 

^ ^'tttatiilt jnyj timei A: aipt^opfi^atlng slg* 

._^ifiB5f|ifii#«^«4:^:^^ 

l*s^i^Qi»5|ii!|B?!^^ on cha 

--,z^zj-i.^ iii^gmiBaisiW'f^^ of . the signal be- 

im.{|if »i^lmi^4i»^<^'>^P^^^ :tha:'time;ctelay.dt>e to 




where w" =^ |wi* w; . . . w,*] and d(w) = [1 e'*"' e'^"^ . . . 
giij-if^Tjn represents the response of the filter to a 
complex sinusoid of frequency u and 6{u) is a vector de- 
scribing the phase of the complex sinusoid at each tap in 
the FIK filter relative to the tap associated with wi. 

Similarly, beamformer response is defined as the ampli- 
tude and phase presented to a complex plane wave as a 
function of location aruJ frequency. Location is in general 
a three dimensional quantity, but often we are only con- 
cerned with one or two dimensional direction of arrival 
(DOA). Throughout the remainder of the paper we do not 
consider range. Fig. 2.2 illustrates the manner in which an 
array of sensors samples a spatially propagating signal. 
Assume that the signal is a complex plane wave with 
□OA e and frequency o). For convenience let the phase 
be zero at the first sensor. This implies xi(k) = e'"* and 
x,(k) = e'*^'^-*"•'^2 < I < J. Ai(e) represents the time delay 
due lo propagation from the first to the I*** sensor. Substi- 
tution into (2.2) results in the beamformer output 



y(k) = e^-" E 2 w,%e - e"*r(e), w) (2.5) 

1-1 p"a 

where — 0. r{$,a}) is the beamformer response and 
can be expressed in vector form as 



r(^,wl = w"d(«,w). 



(2.6) 



The elements of d(^,(u) correspond to the complex 
exponentials e''"^'^'"'**'". In general it can be expressed as 



diB,w) = [1e* 



(2.7) 



where the Ti(fl), 2 :< i < N, are the time delays due to 
propagation and any tap delays from the zero phase refer- 
ence to the point at which the i*^ weight is applied. We 
refer to d(^, w) as the array response vector. It is also 
known as the steering vector or direction vector. Non- 
ideal sensor characteristics can be incorporated into 
d(d, Q>) by multiplying each phase shift by a function 
aAB, ii3>}, which describes the associated sensor response as 
a function of frequency and direction. 
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Thu bf^^nipjitern is defined as the magnitude squared 
of ri0,o}). Nole thai each weight in vv affects both the 
temporal and spatial response of the bcamfcrmer. His- 
toricaily, use of flK filters has been viewed pruvlding 
frequency dependent weights in each channel. This inter- 
prctdtion is accurate but somewhat incomplete since the 
ci>effitients in each filter also influence the spatial filtering 
characteristics n( the beamformer. As a mutti-input single 
output system, the spatial and temporal filtering that oc- 
curs is a result of mutual Interaction between spatial and 
temporal sampling. 

The correspondence between FIR filtering and beam- 
forming is closest when the beamformer operates at a single 
temporal frequency otn and the array geometr>' is linear 
and equi-spaced as illustrated in Fig. 2,3. Letting the sen- 
sor spacing be d, propajjation velocity be c, and U repre> 
sent DOA relative to broadside {perpendicular to the array) 
v\'G have t iO) = (i - 1) (d/c) sin d. In this case wc identify 
the relationship between temporal frequency bi in dM 
(FIR filter) and direction $ in (H0,<*>o) (beamformer) as 
nj - <o^Ui/c.) sin H. Thus, temporal frequency in an FIR fiU 




Figure 2.3 The analogy between an equi-spaced omni- 
dlrBccional narrowband line array and a single-channel FIR 
filter is illustrated in this figure. 



icr corresponds to the sine of direction in a narrowband* 
linear equi-spaced beamformer. Complete interchange of 
beamforming and FIR filtering methods is possible for this 
special case provided the mapping between frequency 
and direction is accounted for. 

The vector notation introduced in (2.6) suggests a vector 
space interpretation of beamforming. This point of view is 
useful both in beamformer design and analysis. We use it 
here in consideration of spatial sampling and array geome- 
try. The weight vector w and the array response vectors 
Of) arc vectors in an N dimensional vector space. The 
angles between w and 6{$.iu) determine the response 
r(B, at). For example, if for some {$, w) the angle between 
w and did, to) is 90^ (i.e., if w is orthogonal to 6{0, to)), then 
the response is zero. If the angle is close to 0'', then the 
response magnitude will be relatively large. The ability to 
discriminate between sources at different locations and/or 
frequencies, say (ffi,aj,) and (0;<,a^?), is determined hy the 
angle between their array response vectors, d{0i <y i) and 
d{0^„oi,} (Cox 119731). 

The general effects of spatial sampling are similar to 
temporal sampling. Spatial aliasing corresponds to an am- 
biguity in source locations. The implication is that sources 
at different locations have the same array response vector, 
e.g., for narrowband sources dfOi^wu) - d(0j,<i>u). This 
can occur if the sensors are spaced loo far apart, tf the 
sensors are too close together, spatial discrimination suf- 
fers as a result of the smaller than necessary aperture; 
array response vectors are not well dispersed in the M 
dimensional vector space. Another t>'pe of ambiguity oc- 
curs with broadband signals when a source at one location 
and frequency cannot be distinguished from a source at a 
different location and frequency, i.e., di<)^,^o^) = d(«3,'tii). 
For example, this occurs in a linear equi-spaced array 
whenever oi, sin 6, = oj^ sin (The addition of temporal 
samples at one sample prevents this particular ambiguity.) 

A primary focus of this paper is on designing response 
via weight selection; however, (2.6) indicates that response 
is also a function of array geometry (and sensor character- 
istics if the ideal omnldirecttonal sensor model is invalid). 
In contrast with single channel filtering where A/O con- 
verters provide a uniform sampling in time, there is no 
compelling reason to space sensors regularly. Sensor loca- 
tions provide additional degrees of freedom in designing 
a desired response and can be selected so that over the 
range of «?, w) of interest the array response vectors are 
unambiguous and well dispersed in the N dimensional 
vector space. Utilization of these degrees of freedom can 
become very complicated due to the multidimensional 
nature of spatial sampling and the nonlinear relationship 
between r(0, cu) and sensor locations. References discus- 
sing array geometry design for response synthesis include 
Unz [19S61, Harrington [19611, Ishimaru [19621, Lo (19631. 
and Skolnik et al. (1964]. 

B. Second Order Statistics 

Evaluation of beamformer performance usually involves 
power or variance, so the second order statistics of the 
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data play an important rofc. We assume the data received 
at the sensors is zero mean throughout the paper. The 
variance or expected power or the beamformer output is 
given by.E{lyj^J = w"E{kx^}w. If the data is wide sense sta- 
tionary, then = E{xx"h <*ata covariancc matrix, is 
independent of time. Although we often encounter non- 
stationary data« the wide sense stationary assumption is 
used in developing statisticalhy optimal beamformers and 
in evaluating steady state performance. 

Suppose X represents samples from a uniformly sam- 
pled lime series having a power spectral density Sioj) and 
no energy outside of the spectral band wbl. Rx can be 
expressed in terms uf the power spectral density of the 
data using the Fourier transform relationship as 

R» » ~" I Sl<i*)d(ft>jd*^'(w)d6> (2,81 

2.TT ^ ' 

with dUo) as defined for {2.4b). Now assume the array data 
X .s due to a source located at direction 9. In like manner 
to the time series case we can obtain the covariance matrix 
of the array data as 

1 

Rx = — I S( widCf?, w)d''( 9, w) 6w . a.9) 

A source is said to be narrowband of frequency w,, ii ran 
be represented as the rank one outer product 

Rv <TZdiH,<oo)d^(B.oto\ (2,10) 

where ai is the source variance or power. 

The conditions under which a source can be considered 
narrowband depend on both the source bandwidth and 
the time over which the source is observed. To illustrate 
this, consider observing an amplitude modulated sinusoid 
or the output of a narrowband filter driven by white noise 
on an oscilloscope. If the signal bandwidth is small relative 
to the center frequency :i.e., if it has a small fractional 
bandwidth;, and the time intevals over whic h the signal is 
observed are short relative to the inverse of the signal 
bandwidth, then each observed wavefonn has the shape 
of a sinusoid. Note that as the observation time interval is 
increased, the bandwidth must decrease for the signal to 
remain .sinusoidal in appearance. It turns out, based on 
statistical arguments, that the observation time bandwidth 
product (IBVVP) is the fundamental parameter that deter- 
mines whether a source can be viewed as narrowband 
sBucklcy 1 19871, Compton (19881). 

An array provides an effective temporai aperture over 
which a source is observed. Fig. 2.2 illustrates this tempo- 
ral aperture 7(0} for a source arriving from direction ft. 
Clearly the TBWP is dependent on the source DOA. An 
array is considered narrowband if the observation TBWP is 
much less than one for all possible source directions. 

Narrowband beamforming is conceptually simpler than 
broadband since one can ignore the temporal frequency 
variable. This fact, coupled with interest in temporal fre- 
quency analysis for some applications, has motivated 
implementation of broadband beamformers with a nar- 
rowband decomposition structure, as illustrated in 



Fig. 2,4, The narrowband decomposition is often per- 
formed by taking a discrete Fourier transform (DFT) of 
the data in each sensor channel using an FFT algorithm. 
The data across the array at each frequency of interest are 
processed by their own beamformer. This is often termed 
frequency domain beamforming. The frequency domain 
beamformer outputs can be made equivalent to the DFT of 
the broadband beamformer output depicted in Fig. 2.1b 
with proper selection of beamformer weights and careful 
data partitioning. This equivalence corresponds to imple- 
menting FIR filters via circular convolution with the DFT. 
C. Beamformer Ciassification 

Beamformers can be classified as either data indepen- 
dent or statistically optimum, depending on how the 
weights are chosen. The weights in a data independent 
beamformer do not depend on the array data and are 
chosen to present a specified response for all signal/ 
interference scenarios. The weights in a statistically opti- 
mum beamformer are chosen based on the statistics of 
the array data to "optimize" the array response. In gen- 
eral, the statistically optimum beamformer places nulls in 
the directions of interfering sources in an attempt to maxi- 
mize the signal to noise ratio at the beamformer output. A 
comparison between data independent and statistically 
optimum beamformers is illustrated in Fig. 2.5. 

The next four sections cover data independent, statisti- 
cally optimum, adaptive, and partially adaptive beamform- 
ing. Data independent beamformer design techniques are 
often used in statistically optimum beamforming (e.g.. 
constraint design in linearly constrained mtnimum vari- 
ance beamforming). The statistics of the array data are not 
usually known and may change over time so adaptive algo- 
rithms are typically employed to determine the weights. 
The adaptive algorithm is designed so the beamformer 



IDKf 



Figur>e S.4 BeajTif orrning is sometimes per*formed in the 
frequency dornain when broadband signals are of Interest. 
This figure iJtustraces trarisf ormatton of the data at each 
sensor into the frequency domain. Weighted combinations 
of the data at each frequency Cbin) are formed. An inverse 
discrete Fourier transform produces the output time 
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Figure 3.1 Part (a) depicts the array configuration and 
defiiies azimuth angle ^ and elevation angle B for a 16 by 
T6 rectangular planar array of sensors. The sensors are 
spaced at one-half wavelength along the x and y axes. The 
narrowband magnitude squared response is illustrated in 
(b) as a function of azimuth and elevation angles. The 
array is phase shift steered to 4> = 30* and 8 = 45'. The 
weights applied to each sensor channel are the product of 
X axis and y axis Dolph-Chebyshev taper weights. 



response converges to a slatislically optimum sulution. 
PartiaDy adaptive beamformers reduce the adaptive algo- 
rithm computational load at the expense of a loss (designed 
to be snialli in .statistical optimality. 



III. DATA INDEPENDENT BEAMFORMINC 

The weights in a data independent beamformer are de- 
signed so the beamformer response approximates a de- 
sired response independent of the array data or data statis- 
tics. This design objective — approximating a desired 
response — is ihe same as that for classical FIR filter design 
(see. for example. Parks and Burrus (1987|). We shall ex- 
ploH iho analogies between beamforming and FIR filtering 
where possible in developing an understanding of the do- 
sign problem and in presenting design procedures. We 
also discuss design problems specific to beamforming. 

Thc^ first part of this section discusses forming beams in 
a classical sense, i.e., approximating a desired response of 
unity at a point of direction and zero elsewhere. Methods 
for designing beamformers having more general forms of 
desired response are presented in the second part. 

i4. Classic ai Beamforming 

Consider the problem of separating a single complex 
frequency component from other frequency components 
using the J tap I IR filter illustrated In Fig. 2.3. If frequency 
to,, is of interest, then the desired frequency response is 
unity at <u« and zero elsewhere. A common solution to this 
problem is to choose w as the vector d(aioi. This choice can 
be shown to be optimal in terms of minimizing the 
squared error between the actual response and desired 
response. The actual response is characterized by a main 
lobe {or beam) and many sidelobes. Since w - diw^,), each 
clement of w ha.s unit magnitude. Tapering or windowing 
the amplitudes of the elements of w permits trading of 
main lobe or 6e.im width against sidelobe levels to form 
the rospor^se into a desired shape. Let T be a I by J diago- 
nal matrix with Ihe real-valued taper weights as diagonal 
elcmenis. The tapered FIR filter weight vector is given by 
Jdiio). A detailed comparison of a targe number of taper- 
ing functions is given in Harris (1978 J, 

tn spatial filtering one is often interested in receiving a 
signal arriving from a known location point B^,. Assuming 
the signal is narrowband (frequency 4y„), a common choi- 
ce for the beamformer weight x'eclor is the array reponse 
vector d{H..,(oJ. The resulting array and beamformer is 
termed a phased array since the output of each sensor is 
phase shifted prior to summation. Fig. 1.2 depicts the 
magnitude of the actual ro.'iponse when w = d(0u,a>i]). As 
tn the f IK filter discussed above, beam width and sidelobe 
levels are the important characteristics of the response. 
Amplitude tapering can be used to control the shape of 
the response, i.e., to form the beam. Figs. 2.5a and b illus> 
trate the effect of amplitude tapering on the response. 

The equivalence of the narrowband linear equi-spaccd 
array and FIR filter (see Fig. 2.3) Implies that the same 
techniques for choosing taper functions are applicable to 



either problem (Schelkunoff 119431). Melhods for choos- 
ing tapering weights also exist for more general array con- 
figurations. If the array is narrowband and the sensors lie 
on a Vine, then methods evolving from continuous spatial 
aperture design can be employed. A desirable smooth 
amplitude distribution function from a continuous aper- 
ture is approximated by a step distribution in a discrete 
aperture (see Ishimaru 119621 for relationship between 
continuous and discrete apertures). 

If the array is planar and factorable, then line array tech- 
niques can be used to synthesize the overall response as 
The product oi two linear array responses. A planar array in 
the xy plane is factorable if its response can be "factored" 
into the product of responses due to line arrays in the x 
and y directions. Fig. 3.1 depicts the beampattern for a 
16 x tfj narrowband planar array. The sensors are spaced 
by one-half wavelength in both the x and y directions and 
the direction of interest is 30*^ bearing and 45" elevation. 
The response is synthesized as the product of Dolph- 
Chebyshcv tapered line arrays in the x and y directions. 

If the beamformer is broadband and employs FIR filters, 
then tapering can be applied independently to both the 
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Figure 3.2 This figure illustrates a dasisical broadband 
beamfarmer. In dessicsr beantf orming (in particular for 
line arrays)! taper waiglits {p^*.t^r"* tj} are used to 
shape the bQampattam. To stieer Q bpoodband enrey. fil- 
ters are used co compensace for isropagatlon delays of 
the desired source. 



sensor outputs and FIR fillers as illustrated in Fig. 3.2. The 
lapcr weights are chosen to shape the spatial response 
and the FIR filter coefficients to present a desired temporal 
response. As noted in Section II. spatial and temporal re- 
sponse interact so the spatial and temporal responses 
cannot be synthesized completely Independently. One ex- 
ample of where the structure of Fig. 3.2 is used is delay 
sum beamforming. Here the FIR filters approximate the 
propaj^alion delays (linear phase over the frequency band 
of interest) and the taper weights are chosen to shape the 
main beam and sidetobe structure of the spatial response, 
ff. Genera! Data fndepenthnt Respqase Design 

The methods discussed in this subsection apply to de- 
sign of bcamformers that approximate an arbitrary desired 
response. This is of interest in several different applita- 
tions. For example, we may wish to receive any signal 
arriving from a range of directions, in which case the de- 
sired response is unity over the entire range. As another 
example, we may know that there is a strong source of 
interlcrence arriving from a certain range of directions, in 
which case the desired response is zero in this range. 
These two examples are analogous to bandpass and band* 
stop FIR filtering. Although we are no longer "forming 
beams/' it is conventional to refer to this type of spatial 
filter as a beamformer. 

Consider choosing w so the acluat response r(tf. w) ~ 
w"d(fl, w) approximates a desired response r,j(!>, w). 
Ad hoc techniques similar to those employed in FIR filter 
design can he used for selecting w; however, here we 
only constder choosing w to minimize the weighted Lp 
norm of the difference between desired and actual re- 
sponse. Weighted Lp approximation is utilized in several 
establi.'ihfxj FIR filler design techniques. The most com- 
monly used norms are (minmax) and I2 (least squares). 
Specific techniques include (see Parks and Burrus 119871): 

1) Windowing of an ideal filter's unit pulse response 
(minimizes L2 norm over continuous «); 

2) Frequency response sampling and linear weighted 
least squares (minimized L2 norm over discrete w); 



3) Minmax design with the Remez exchange algorithm 
(minimizes L. norm over discrete w); 

4) Minmax complex and magnitude response design 
(minimizes L« norm over discrete oj). 

FIR filter design corresponds to a polynomial approxi- 
mation problem since the frequency response (2.4b) is the 
discrete Fourier transform of the FIR filter weight se- 
quence. Several of the above methods exploit this poly- 
nomial structure. 

Excluding the cases for which beamformer design can 
be reduced to equi-spaced line array geometries, beam* 
former design is not a polynomial approximation problem. 
In general, the response in (2.6) is a weighted sum of 
exponentials raised to non-integer powers. Thus, the L» 
methods (3 and 4) arc not applicable since they are based 
on the alteration theorem of polynomial approximation. 
The windowing method (1) is based on the discrete time 
Fourier transform and is also not applicable. However, the 
L2 procedure using linear weighted least squares (2) is 
applicable. 

To illustrate data independent beamformer design via 
optimization, consider minimizing the squared error be- 
tween the actual and desired response at P points ufi), 
1 i P. If P > N, then we obtain the overdctermined 
least squares problem 

min jA"w - T^\^ (3.1) 

where 

A = [dlDi.oj,) dirtz, III?) -deep, wp)j; 
r,i = [r«r(ei,w,) rd(fl2,w2)"-rd«?p, wp)r. 

Provided AA" h invertible (i.e., A is full -rank), then the 
solution to (3.1) is given as 

w = A'fd (3.2) 

where A* = (AA^)" *A is the pseudo inverse of A. Fig, 3.3 
depicts the response of a beamformer design using (3.2). 

A note of caution is in order at this point. The white 
noise gain of a beamformer is defined as the output power 
due to unit variance white noise at the sensors. Thus, the 
norm squared of the weight vector, w*V, represents the 
white noise gain. If the white noise gain is large, then 
the accuracy by which w approximates the desired re- 
sponse ts a moot point, since the beamformer output will 
have a poor SNR due to white noise contributions. If A is 
ill-rondittoned, then w can have a very large norm and still 
approximate the desired response. The matrix A is ill- 
conditioned when the numerical dimension of the space 
spanned by the dii)\<ol. 1 i P, is less than N. For ex- 
ample, if only one source direction is sampled, then the 
numerical rank of A is approximately given by the TBWP 
for that direction. Low rank approximates of A and A* 
should be used whenever the numerical rank is less 
than N. This ensures that the norm of w will not be unnec- 
essarily large. 

Specific directions and frequencies can be emphasized 
in (3.1) by selection of the sample points Wi.tod and/or 
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Figure 3.3 This plat da plctis a independent beam- 
pattern evaluated ac eig^it f pequencies tin the nprmsliied 
frequency interval [2ir/5, 47r/53:- Therweiahts are de- 
signed according co (3; 'VJ ^Arfth adpsiredresponseof tmrty 
gain and tinear phase- ovar ISW/S. 47r/5J at a djrfection vf 
16 degrees. The array is lUnaar etM-^P^9^ y'li^ .sixtee 
sensors spaced at. ona-'half w^l«ngiii tfqr frequency 
4v/5 and Five tap RR filjter8;'a(^.ijlBad bi-each aensor 
channel. 7 • 



unequally weighting of the error at each {Oi, wi). Parks and 
Burrub 119871 discuss this in the context of FIR filtering. 
Kumar and Murthy 11977J consider unequal weighting of 
the error lo obtain minimax response design for beam- 
former wei^^hts rn a linear narrowband array. In general, 
guidelines lor selection of the error weighting and to.) 
are not available. 



There are several alternatives lo Lp optimization for 
general data independent response design, including 
methods discussed in Boiler and Unz (1967J and Sarkzgiri 
and Butler [1971]. 



IV. STATISTiCAUY OPTIMUM BEAMFORMINC 

In statistically optimum beamforming the weifihis are 
chosen based on the statistics of the data received at the 
array. Loosely speaking, the goal is to "optimize" the 
beamformcf response so the output contains minimal 
contributions due to noise and signals arriving from direc- 
tions other than the desired signal direction. We discuss 
below several different criteria for choosing statistically 
optimum beamformer weights. Table 4.1 summarizes 
these different approaches. Where possible, equations 
describing the criteria and weights are confined to 
Table 4.1 . Throughout the section we assume that the data 
is wide-sense stationary and that its second order statistics 
are known. Determination of weights when the data statis- 
tics are unknown or time varying is discussed In the fol- 
tov^nng section on adaptive algorithms. 

A. Multiple SiMobe Canct^tler 

The multiple sidclobe canceller (MSC) is perhaps the 
earliest statistically optimum beamformer. An MSC con- 
sists of a "main channel" and one or more *auxillary 
channels" as depicted in Fig. 4.1a. The main channel can 
be either a single high gain antenna or a data independent 



TABLE 4.1 
' ::: i Summary: OP. Optimum BEAMFonMERS 



Type 
Definitions 



Criterion 

Optimum 
Weights 

Advantages 



Disadvantages 



...:.M5C.:. . 
xi ^ auxi Itary; iciata ; . . 
ym ^primary, dita 

outptit: Yf'y^w^^i 

min E{|y„-w?K«P} • 

Simple 



Requires absence of 
desired signal from 
auxiliary channels for 
weight determiriation 



References Applebaum 11976) 



jieference Signal 
X— array data 
y^- — desired signal 
rHd^E{«y3} 

output:. y=w"x 



mjn E{|y-yrfi='} 

Direction of 
desired signal 
can be unknown 

Must generate 
reference signal 



WIdrow (1967] 



MaxSNR 

x«s+n^ — array data 
s — signal comporient 
h — noise component 
R.=E<s$"} 
R„-E{nn"} 
output: y=w'^x 

w"R.w 

max u 

W W RnW 
Rn 'R»W— A W 

True maximization 
of SNR 

Must know R» 
and R„, 

Solve generalized 
eigenproblem for 
weights 

Monzingo and 
Miller [19801 



LCMV 
X — array data 
C — constraint matrix 
i — response vector 
Rx-E{xx"} 
output: y=w"x 



minf w"R,w}s.t .C^w = f 
w 

w-R.-X[c^Rr'c)-'f 



Flexible and general 
constraints 



Computation of 
constrained 
weight vector 



Frost [19721 
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Figure 4.1 The multiple sidelobe canceller (MSCl con> 
sists of a main channel and several auxiliary channels as 
illuscrated in a). The auxiliary channel weights are chosen 
to "cancel" interference entering through the etdelobes of 
the main channel response. Part W depicts the main chan- 
nel, auxiliary branch, and overali system response when an 
interferer arrives from direction Gi. 



bcamtormcr (see Section III). It has a highly directional 
response, which is pointed in the desired signal direction. 
Interfering signals arc assumed to enter through the main 
channel ssdclobes. The auxiliary channels also receive the 
interfering signals. The goal is to choose tfie auxiliary 
channel weights to cancel the main channel interference 
component. This implies that the responses to interfcrers 
of the main channel and linear combination of auxiliary 
channels must be identical. The overall system then has a 
response of zero as illustrated In Fig. 4.1b. In general, 
requiring ^ero response to all interfering signals is either 
not possible or can result In significant white noise gain. 
Thus, the weights are usually chosen to trade off inter- 
ference suppres.sion for white noise gain by minimizing 
the expected value of the total output power as Indicated 
in Table 4.1. 

Choosing The weights to minimize output power can 
cause cancellation of the desired signai, since it also con- 
tributes to total output power, in fact, as the desired signal 
gets slronger it contributes to a larger fraction of the total 
output power and the percentage cancellation increa.ses. 
Clearly ihis is an undesirable effect. TheMSC is very effec- 
tive in applications where Ihe desired signal is very weak 
(relative to the inlerference), since the optimum weights 
will not pay any attention to it, or when the desired signal 
is known to be absent during certain lime periods. The 
weif>hts can be adapted in the absence of the desired 
signal and frozen when it is present. 
B. Use of Reference Signa! 

If the desired signal wore known, then Ihe weights could 
be <:ho&en to minimize the error between the beamformer 
output and Ihe desired signal. Of course, knowledge of 
the desired signal eliminates the need for beamforming. 
However, for some applications enough may be known 
about the desired signal to generate a signal that closely 



represents it. This signal is called a reference signal. As 
indicated in Table 4.1, the weights are chosen to minimize 
the mean square error between the beamformer output 
and the reference signal. 

The weight vector depends on the cross covariance be- 
tween the unknown desired signal present in x and the 
reference signal. Acceptable performance is obtained pro- 
vided this approximates the covariance of the unknown 
desired signal with Itself. For example. If the desired signal 
is amplitude modulated, then acceptable performance is 
often obtained by setting the reference signal equal to 
the carrier. It is also assumed that the reference signal is 
uncorrrelated with interfering signals in x. The fact that 
the direction of the desired signal does not need to be 
known is a distinguishing feature of the reference signal 
approach. 

C. Maximaation of Signal to Notse Ratio 

Here the weights are chosen to directly maximize the 
signal to noise ratio (SNR) as indicated In Table 4.1. A gen- 
eral solution for the weights requires knowledge of both 
the desired signal, and noise, R„, covariance matrices. 
The aftainabilily of this knowledge depends on the appli- 
cation. For example, in an active radar system Rn can be 
estimated during the time that no signal is being trans- 
mitted and Rs can be obtained from knowledge of the 
transmitted pulse and direction of interest. If the signal 
component is narrowband, of frequency o/, and direction 
e, then = cr'^diti, uf)d"{0, di) irom the results in Sec- 
tion IL In this case the weights are obtained as 

w = aR„'d(0,w) <4.1) 

where the a Is some non-zero complex constant. Substi- 
tution of (4.1) Into the SMR expression shows that the SNR 
is Independent of the value chosen for u. 

D. Linearly Constraint Minimum Variance Beamforming 

in many applications none of the above approaches is 
satisfactory. The desired signal may be of unknown 
strength and may always be present, resulting in signal 
cancellation with the MSC and preventing estimation of 
signal and noise covariance matrices In the maximum SNR 
processor. Lack of knowledge about the desired signal 
may prevent utilization of the reference signal approach. 
These limitations can be overcome through the applica- 
tion of linear constraints to the weight vector. Use of lin- 
ear constraints is a very general approach that permits 
extensive control over the adapted response of the 
beamformer In this subsection we illustrate how linear 
constraints can be employed to control beamformer re- 
sponse, discuss the optimum linearly constrained beam- 
forming problem, and present the generalized sidelobe 
canceller structure. 

The basic idea behind linearly constrained minimum 
variance (LCMVl beamforming Is to constrain the re- 
sponse of the beamformer so signals from the direction of 
interest are passed with specified gain and phase. The 
weights are chosen to minimize output variance or power 
subject to the response constraint. This has the effect of 
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preserving the desired signal while minimizing con- 
tributions to the output due to interfering signals and 
noise arriving from directions other than the direction of 
interest. The analogous FIR filter has the weights chosen 
to minimize the filler output power subject to the con- 
straint that the filter response to signals of frequency <uo 
be unity. 

In Section II we saw that the beamformer response to a 
source at angle 6 and temporal frequency ut is given by 
w'*d«l..a/). Thus, by linearly constraining the weights to 
satisfy w"d(0,w) = g, where g is a complex constant, wc 
ensure that any signal from angle 0 and frequency cj is 
passed to the output with response g. Minimization of 
contributions to the output from Interference (signals not 
arriving from $ with frequency o)) is accomplished by 
choosing the weights to minimize the expected value of 
output power or variance E{|y|"} - w"R.w. The LCMV 
problem for choosing the weights is thus written 

'f w*'R^w subject to d'^«?,w)w = g*. (4.2) 

The method ot Lagrange multipliers can be used to solve 
(4.2) resulting in 

Note that in practice the presence of uncorrelatcd notsc 
will ensure that is invertiWe. If g - 1, then (4.H) is often 
termed the minimum variance distortionless response 
(MVDR) beamformer. It can be shown that (4.3) is equiva- 
lent to the maximum SNR solution given in (4.1) by substi- 
tut'm\» (r^d{e,(o)&\o.<iA + Rr. for R, cn (43) and applying 
the matrix inversion lemma. 

The single linear constraint in (4.2) is easily generalized 
lo multiple linear constraints for added control over the 
beampattern. For example, if there is a fixed interference 
source at a l<nown direction th, then it may be desirable to 
forte /.ero gain in that direction in addition to maintaining 
the response g to the desired signal. This is expressed as 

K::;]"-[v]- 

If there are L < N linc^ir constraints on we write them 
in the form C"w = f w})ere the N by L matrix C and L 
dimensional vector f arc termed the constraint matrix and 
responsR vector. The constraints are assumed to be lin- 
early independent so C has rank L. The LCMV problem and 
soUition with this more general constraint equation are 
given in Table 4.1 . 

Constrainl Design. Several different philosophies can 
t)e employed for choosing the constraint matrix and re- 
sponse vector. Wo discuss point (Kelly and Lovin (1%4|). 
derivative (Qwsloy [19731. Er and Cantoni 11983]), and ei- 
genvector (Buckley |1987]) constraints below. In many 
applications, a combination of the different types of con- 
straints is most effective. Each linear constraint uses one 
degree of freedom in the weight vector so with L constraints 
there are only N • L degrees of freedom available for 
minimising variance. 



Point constraints fix the beamformer response at points 
of spatial direction and temporal frequency. Equation (4.4) 
represents an example of two point constraints on w. The 
number of points at which response can be constrained is 
limited to N. If N constraints are used then there are no 
degrees of freedom left for power minimization and a data 
independent beamformer is obtained. 

Derivative constraints arc employed to influence re- 
sponse over a region of direction and/or frequency by 
forcing the derivatives of the beamformer response at 
.some point of direction and frequency to be zero. They 
are usually employed in conjunction with point con- 
straints. An example where derivative constraints are use- 
ful is when the desired signal direction is only known 
approximately. If the signal arrives near Ihe direction at 
which a point constraint is employed, then application of 
a derivative con.straint at that point prevents the beam- 
former from synthesizing a response off zero lo the de- 
sired signal. 

Eigenvector constraints are based on a least squares ap- 
proximation to the desired response and are typically used 
to control beamformer response over regions of direction 
and/or frequency. Constraining the beamformer response 
in a least squares sense ensures that the mean square error 
between desired and actual beamformer response over a 
region is minimized for a given number of constraints. In 
this sense eigenvector constraints are efficient. Consider 
designing a set of constraints which will control the beam- 
former response to a source from direction 6^ over the 
frequency band wt,]. The dimension of the span of 
6iO<„ to) over this band is approximately given by the 
source TBWP discussed in Section II. Eigenvector con- 
straints are derived from (3.1) by choosing P significantly 
greater than the TBWP. The faj» then oversampic lai„6ih| 
and A is ill-conditioned. A rank L approximation of A is 
obtained frotn its singular value decomposition 

Al==V3:^U" (4.5) 

where Xi is an L by L diagonal matrix containing the largest 
singular values of A, and the L columns of V and U are 
respectively the left and right singular vectors of A corre- 
sponding to these singular values. Replacing A in (1.1) by 
its rank L approximate (4.5) and bringing US, to the right 
side (the pseudo inverse of U is U'*), yields 

V"w^iVU"rd. (4,6) 

Equation (4.6) has the same form as the constraint equa- 
tion C"w *= f. The columns of V correspond to the eigen- 
vectors of AA"; hence the name eigenvector constraints. 
(Note that AA^ represents an approximation of in (2.9) 
itS(w) = 1.) 

fxampfe. Fig. 4.2 depicts the response of an LCMV 
beamformer at eight equally spaced frequencies on the 
interval f2tr/5,4ir/5| when two interferers arrive from 
-17.5 and —5.75 degrees in the presence of white noise. 
The interferer power levels relative to the white noise are 
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-liB -05 ■ -ft} • -«>• HZO' 0' ap <0 (O Vf iX^ 

Figure 4;2 This plot deprcte the LCMV beamformer 
magnitucle riasppnse.at eigtit frequencies on the normah 
ized fpsquehcY=int9pyBt t2WS.4<ir/5J whan two intsp- 
ferers arrive from directions 5. 76" and - 1 7.5 degrees 
in t*T8: presence of white ndse. The mterfarera have a 
white spectrum on t2Tr/5.4tr/5] and have powers of 40 
and" 30 dB relative. KD the White noise, respectively. The 
constraihts are ide^cpied according, to (4.G) to present 
unit gain and linear phae^ over (Str/S^ 49r/5J at a OOA of 
IB defin^es. The array je. the same as that of Fig. 3.3. 



decomposition can he used to represent any w. Since 
C"y - 0, we must have 



30 and 40 dB, respectively, and the inlerferers havu a flat 
spectrum on l2;r/5,47r/S]. The array ha.s sixteen linear 
equi-spaccd sensors with five taps per sensor. The tap 
spacing Is normalized to one second and the sensors are 
spaced at one-halt wavelength corresponding to fre- 
quency 4-^/5. The response is constrained to pass signals 
arriving from 18 degrees in the band [2ir/5,47r/5J, with 
unit gain and linear phase using ten ergcnvector con- 
straints designed from (A.b). The effectiveness of the con- 
straints is evident, since all the frequency curves pass 
through zero dH at 18 degrees. The response has nulls in 
the directions of the interferers with the deeper null corre- 
sponding to the stronger interferer. The response as a 
function of frequency for the interferer directions is plot- 
ted in f ig. 6.2. The array gain is SO dti for this example. 

Genera/izecf S'tdehbe Canceffer. The generalized slde- 
lobc canceller tCSO represents an alternative formulation 
of the I.CMV problem^ which provides tn&ight, is useful for 
analysis, and can simplify LCMV beamformer implementa- 
tion. It also illustrates the relationship between the MSG 
and ICMV boamforming. Essentially, the GSC is a mecha- 
nism for changing a constrained minimization problem 
into unconstrained form. Perhaps the first reference to 
this concept is in Hanson and Lawson [1%9J, where a pro- 
cedure for transforming constrained least squares prob- 
lems U) uncontitrained lea.st ."iquarcs problems is given. 
Griffiths and Mm [1982] applied the same concept to LCMV 
beamforming and coined the term CSC Similar tech- 
niques were discussed in Applcbaum and Chapman (1976] . 

Suppose wc decompose the weight vector w into hvo 
orthogonal components w„ and -v (w = w,, - v) thai lie 
in the range and null space of C. respectively. The range 
and null space of a matrix span the entire space so this 



Wo « ClC'*C)-'f 



(4.7) 



if w is to satisfy the constraints. (4.7) is the minimum L> 
norm solution tothe underdetcrmined equivalent of (3.1). 
The vector v is a linear combination of the columns of an 
N by N-L matrix Cr,(v = Cw^) provided the columns of C„ 
form a basis for the null space of C. C„ can be obtained 
from C using any of several orthogonalization procedures 
such as Cram-Schmidt, QR decomposition, or singular 
value decomposition. The weight vector w = w„ - C^Wn 
is depicted in block diagram form in Fig. 4.3. The choice 
for and C> implies that w satisfies the constraints inde- 
pendent of Wn and reduces the LCiVlv problem to the 
unconstrained problem 



The solution is 



CnWn]**RB[Wo d.WnJ. 



w.. - (C|^R.C„)-'C:.'R.w„ 



{4.6) 



(4.<J) 




Figure 4.3 The generalized sidelobe canceller CGSC) rep- 
rasents an tmplementation of the LCfVIV beamformer in 
which the adapcivB weights are unconstrained. It consists 
of a fixed tteamfonmeri Wo, signal bicx:kfng matrix. Cm. and 
adaptive weight vector «f„. 



The primary implementation advantages of this alter- 
nate but equivalent formulation stem from the facts that 
the weights Wn are unconstrained and a data independent 
beamformer w„ is implemented as an integral part of the 
adaptive beamformer. The unconstrained nature of the 
adaptive weights permits much simpler adaptive algo- 
rithms to be employed and the data independent beam- 
former is useful in situations where adaptive signal 
cancellation occurs (see subsection E). 

As an example, assume the con.straints are as given in 
f4.2). (4.7) implies w., - g*d(^. fc;)/|d"(f>. (u)d(#, w)]. C„ sat- 
isfies d*\ti,ii}}Cn - 0 so each column IC^).; 1 <: i < N'^L. 
can bo viewed as a data Independent beamformer with a 
null in direction 6 at frequency o) \ d"(0, w) ICL = 0. Thus, 
a signal of frequency u and direction 6 arriving at the ar- 
ray wilf be blocked or nulled by the matrix C. In general, 
if the constraints are designed to present a specified re- 
sponse to signals from a set of directions and frequencies, 
then (he columns of Cn will block those direc:tions and 
frequencies. This characteristic has led to the term "block- ' 
ing matrix" for Cn- These signals are only processed by w^ 
and since w« satisfies the constraints, they are presented 
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with the desired response independent of w„. Signals 
Ifom directions and frequencies over which the response 
is not constrained will pass through the upper branch in 
Fr^. 4.3 with some response determined by w^. The lower 
branch chooses Wn to estimate the signals at the output of 
Wo as a linear connbination of the data at the output of the 
blocking matrix. This is similar to the operation of the MSC. 
in which wiiiyhts are applied to the output of auxjilary 
sensors in order to estimate the primary channel output 
{see Fig. 4.1). 

F. Si'xnaf Canceltadon in StatisticaUy Opiimum 

Beamforming. 

Optimum beamiorming requires some knowledge of 
the desired signal characteristics, either its statistics (for 
maximum SNK or reference signal methods), its direction 
(lor the MSG), or its response vector d(tf, u/) (for the LCMV 
bcamformer). If the required knowledge is inaccurate^ the 
optimum beamformer will attenuate the desired signal as 
if it were interference. Cancellation of the desired signal is 
often significant, especially if the SNR of the desired signal 
is large (C^ox [197.*^]). Several approaches have been sug- 
gested to reduce this degradation (e.g., labion 11986], Cox 
et al. : 19871). 

A second cause of signal cancellation is correlation be- 
tween the desired signal and one or more interference 
signals. This can result either from muUipath propagation 
of a desired signal or from smart (correlated) Jamming. 
When interference and desired signals are uncorrelated 
the beamformer attenuates interferers to minimize output 
power. However, with a correlated i nterf erer the beam- 
former minimizes output power by processing the inter- 
fering signal in such a way as to cancel the desired signal. 
If the inlerferer is partially correlated with the desired 
signal, then the beamformer will cancel the portion of the 
desirtfd signal that is correlated with the interferer. Meth- 
ods for reducing signal cancellation due to correlated 
interference have been suggested (e.g., VVidrow ct al. 
[1982J, Shan and Kailath (1985], Yang and Kaveh 11987]). 

V. ADAPTIVE ALGORITHMS FOR BEAMFORMINC 

The optimum beamformer weight vector equations 
listed in Table 4.1 require knowledge of .second order sta- 
li.stics. These statistics are usually not known, but with the 
assumption of ergodlcity, they (and therefore the opti- 
mum weights) can be estimated from available data. Sta- 
tistics may also change over time (e.g., due to moving 
interferers). To solve these problems, weights are typically 
determined by adaptive algorithms. 

There are Uvo basic adaptive approaches: 11 block adap- 
tation, where statistics are estimated from a temporal 
block of array data and used in an optimum weight equa- 
tion; arKl 2) continuous adaptation, where the weights are 
adjusted as the data is sampled such that the resulting 
weight vector sequence converges to the optimum solu- 
tion. If a nonstalionary environment is anticipated, block 
adaptation can be used, provided that the weights are 
recomputed periodically. Adams, et al. I1980J and others 



have described applications of block data processing. 
Continuous adaptation is usually preferred when statistics 
are time-varying or (for computational reasons) when the 
number of adaptive weights ^f is moderate to large (values 
of M > 50 are not uncommon). 

Among notable adaptive algorithms proposed for 
beamforming are the Howells-Applebaum adaptive loop 
developed in the late 195U's and reported by Howells 
[1*i66,1976] and Applebaum 1196b], and the Frost LCMV 
algorithm (19721. Rather than recapitulating adaptive algo- 
rithms for each optimum beamformer listed in Table 4.1 
(for this see texts by Monzingo and Miller (1980], Hudson 
119811 and Compton (1988J), we take a unifying approach 
using the standard adaptive filter problem illustrated in 
Fig. 5.1. 




Figure 5.1 The standard adaptive filter configuration. 



In Fig. 5.1 the weights are chosen to estimate the de- 
sired signal Yd as linear combination of the elements of the 
data vector u. We select to minimize the MSE 

trit - v/^r^f ~ T^Wm -i- WmRwWm , (5.1) 

where ai = f {jy^j^}. r„„ « E\uyrA and R., - E{uu"]. iS.l) is 
minimized by 

w,,„ = R-V„tf. (5.2) 

Comparison of (5.1) and the criteria listed in Table 4.1 indi- 
cates that this standard adaptive filter problem is equiva- 
lent to both the MSC beamformer problem (with y</ = y^ 
and u ^ X;,) and the reference signal beamformer problem 
(with u = x). The LCMV problem is apparently different. 
However closer examination of Figs. 4.3, 5.1, and Equa- 
tions (4.9), (5-2) reveals that the standard adaptive filter 
problem is equivalent 10 the LCMV problem implemented 
with the CSC structure. Setting u « C»x and = wi,'x 
implies Rp = CVRkC„, and r^a = Cj^R.Wo- The maximum 
SNR beamformer ran not In general be represented by 
Fig. 5.1 and Equation (5.2). However, it was noted 
after (4.3) that if the desired signal is narrowband, then 
the maximum SNR and the LCMV beamformers are 
equivalent. 

The block adaptation approach solves (5.2) using esti- 
mates of R„ and r^^ formed from K samples of u and y<, : 
lilk). yAk); 0 ^ k ^ K - ^. The most common arc the 
sample covariance matrix 

' ^ 1 * ' 

77 S u(/c)u"(*J 
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and sample cross-covariance vector 

Pertormance analysis and guidelines for selecting ihe 
block size K are provided In Reed et al. 119741. 

Continuous adaptation afgorithms arc easily developed 
in terms o! Fig. 5.1 and Equation {5.1 >. Note that /(wm) 
is a quadratic error surface. Since the quadratic surface's 
"Hessian'' R„ is the covariancc matrix of noisy data, it is 
positive definite. This implies that the error surface is a 
"bowl." The shape of the bowl is dcicrmined by the cigen* 
structure or R„. The optirDum weight Wop* corresponds to 
the bortom of the bowJ. 

One approach to adaptive filtering is to envision a point 
on the error surface that corresponds to the present 
weight vector vcxi\k). We select a new weight vector 
v/Uk + 1) so as to descend on the error surface. The gra- 
dient vector 

« .1 

(JWai ,v*m ****** I 

- -2r^ + 2R.;W«(/c» (5.3) 

tells us the direction in which to adjust the weight vector. 
Steepest descent, i.e., adjustment in the negative gradient 
direction, leads to the popular least-mean-square (LMS) 



adaptive algorithm. The LMS algorithm replaces V„^t, 
with the instantaneous gradient estimate 
-2\M{k)Y:,{k) - u(^)u"(Jt)Ww(fc>I. Denoting y{k) « y^/M) " 
v^sMuik). we have 



(5.4) 



The gain constant ^ controls convergence characteristics 
of the random vector sequence w^^f^f). Table 5.1 provides 
guidelines tor its selection. 

The primary virtue of the LMS algorithm is its simplicity. 
Us performance i.s acceptable in many applications; how- 
ever, its convergence characteristics depend on the shape 
of the error surface and therefore the eigen structure of R„. 
When the eigenvalues arc widely spread, convergence 
can be slow and other adaptive algorithms with better 
convergence characteristics should be considered. Alter- 
native procedures for searching the error .surface have 
been proposed in addition to algorithms based on least- 
squares and Kalman filtering. Roughly speaking, these 
algorithms trade-off computational requirements with 
speed of convergence to w,^,. We refer you to texts on 
adaptive filtering for detailed descriptions and analysis 
(Widrovv and Stearns 11985), Haykin [1986]. Alexander 
11986], Treichler et at. [1987]. and others). 

One alternative to LMS is the exponentially weighted 
recursive least squares (RLS) algorithm. At the K'^ time 



TABLE 5.1 

COMPAWSON OF THE LMS AND RLS WfefGHT ADAPTATION ALGORITHMS 



Algorithm 
Initialization 



Update 
liquations 



Multiplies 
per update 

Performance 
Characteristics 



IMS 

0. < < ; ^ 



Trace [R.,] 



2M 



Under cerlalnr conditions^ convergence of 
WM(/r):to the statisticajlly optimum Weight' 
vector Wi^: in the mean-square: sense b 
guaranteed Hf fx; is. chosen as indicated 
above, The convergence rate is govetrred 
by the eigenvalue spread ofR«, For Ikijge 
eigenvalue spiread, convergence ca/t be . 
very:s(6W;:-. ■ ' : " ■ . 



RIS 
wv(0) = 0 
P(0) ^ fi-^^l 

6 small r I identity matrix 

v<k) - ?{k - 1)u(Af) 

^^^^ " 1 + A-u'Wv(ik) 
a{k) - yJik) - v^{k - ^)u{k) 
w^X*) = vnt^ik - 1) + V{k)a*ik) 
?{k) = A-V(fc ~ 1) - k-'Vik)^f*'{k) 

AM^ + 4/W + 2 



The YftAk) represents the least- 
squares solution at each instant 
■k and are optimum in a deter- 
ministic sense. Convergence to 
the statistically optimum weight 
vector w«pi is often faster than 
that obtained using the IMS 
algorithm because it is inde- 
pendent of the eigenvalue 
spread of Ru. 
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step, w,JK) Is chosen to minimize a weighted sum of past 
squared errors 

min X * y^M - wZ{K)v{k)\'. (5.5) 

A is a positive constant less than one which doiermines 
how quickly previous data are deemphasized. The RL5 
afgorithm is ohlained from (5.5) by expanding the mag- 
nitude squared and applying the matrix inversion Jemma. 
Tabic 5.1 summarizes both the IMS and RLS algorithms. 

VI. INTERFERENCE CANCELLATION AND PARTIALLY 
ADAPTIVE BEAMFORMINC 

The computational requirements of each update in 
adaptive alKorithm.<s are proportional to either the weight 
vector dimension M or dimension squared <M'). If M is 
large, this requirement Is quite seven: and for practical real 
time implementation 11 is often necessary to reduce M. 

The expression deferens of freedom refers to the num- 
ber of unconstrained or "free" weights In an imple- 
mentation. For example, an LCMV beamformer with L 
constraints on N weights has N-L degrees of freedom; the 
GSC implementation separates these as the unconstrained 
weight vector w„. There are M degrees of freedom in the 
structure of Hig. 5.1. A fully adaptive beamformer uses all 
available degrees of freedom dLn63LpartiaUy adaptive beam- 
former uses a reduced set of degrees of freedom. Re- 
ducing degrees of freedom lowers computational 
requirements and often improves adaptive response 
time,' However, there is a performance penalty associated 
with reducing degrees off freedom. A partially adaptive 
beamformer cannot converge to the same optimum solu- 
tion as the fully adaptive beamformer. The goal of partially 
adaptive beamformer design is to reduce degrees of free- 
dom without significant degradation in performance. 

The discussion in this section Is general, applying to 
different types of beamformors although we borrow much 
of the notation from the CSC. VVc assume the beamformer 
l.sdnscribed by the adaptive structure of Fig. 5.1 where the 
desired signal yj is obtained as ya - w^x and the data 
vector u as u — T"x. Thus, the beamformer output Is 
y w"x where w *- w„ - Twm. In order to distinguish be- 
tween fully and partially adaptive implementations wc 
decompose T into a product of two matrices ChTm. The 
definition of C- depends on the particular beamformer 
and Tm represents the mapping which reduces degrees of 
freedom. The MSC and CSC arc obtained as special cases 
Of this representation. In the MSC Wo is an N vector that 
selects the primary sensor, C^ is an N by N-1 matrix that 
selects the N-l possible auxiliary sensors from the com- 
plete set of N sensors, and Tm is an N-1 by M matrix that 
selects the M auxiliary sensors actually utilized. In terms of 

'Adaptive alj^orithm convergence thnmcieristics have not been dis- 
cussed in ihi<> paper. Generally, more data are required to derive an 
accurate esiimaie of a larger optimum weight vector with biock adap- 
tive processing nr RLS |Keed Gt at., 1974J. With IMS. convergence is 
governed by the covariance matrix eigenvalue spread, which tcndf to 
be larger for larger dimensional problems. 



the CSC, w« and Cn are defined as in Section IV and Tv, is 
an N-L by M matrix that reduces degrees of freedom 

(M < N-U. 

This section begins by considering the interference 
cancellation process In these general beamformer imple- 
mentations. This develops the intuition required for under- 
standing why and how the number of adaptive weights can 
be reduced. We conclude this section by surveying differ- 
ent partially adaptive beamformer design philosophies. 
A, interference Cancelfation Vs Degrees of freedom. 

The results in this subsection depend on T and are Inde- 
pendent of the individual terms C» and Tm. We assume 
that the beamformer dues not cancel the desired signal 
(sec Section IV.E.) and that the optimum weights affect 
only interferers and uncorrclated noise. This .simplifies 
the analysis by permitting us to exclude consideration of 
the desired signal. 

Suppose a narrowband interfering source of frequency 
Wo arrives at the array from dircclion 9,, The response of 
the Wo branch is g, = w^'dCf^i. w J. Perfect cancellation or 
this source requires w"d((?,,<i/,.) = 0 so we must choose 
wm to satisfy 

\v\:j"6i9,,ioJ = g,. (6.1) 

If we asume that T^df^i. Wo) is nonzero. (6.11 represents a 
system of one equation in M unknowns (elements of Wm) 
for which a solution always exists. To simultaneously can- 
cel a second interferer located at <?3, Ww must satisfy 

wr,[T^d«J,, oo) T"d(fl,, - Lgi gil <6.2) 

where gj = w^^Cd^.^vK Assuming T^dCfl^, w„J and 
T^tHiOi. w,J are linearly independent and nonzero, and 
provided M > 2, then at least one Wm exists that satisfies 
(6.2). Continuing this reasoning, we see that Wm can be 
chosen to cancel M narrowband interferers (assuming the 
T"d(0,. Wo) arc linearly independent and nonzero), inde- 
pendent of T. Total cancellation occurs if wm is chosen so 
the response of Twm pcrieclly matches the w„ branch re- 
sponse to the interfc'rers. 

So far we have only considered narrowband point inter- 
ferers. Uncorrelated noise will be present in any real sy.s- 
tern and contributes to the output power. In an optimum 
beamformer Wm is chuscn (o minimize the overall output 
power. Recall that the output power due to uncorrelated 
noise is proportional to the L:^ norm squared of the overall 
weight vector w (white noise gain). The norm of w can 
become large when Wm is chosen to provide total inter- 
ference cancellation, depending on the choice for T and 
the interferer locations. Thus, although in principle point 
sources of energy in direction and frequency can be totally 
canceled with one weight per interferer Independent of T, 
the presence of uncorrelated noise results in the degree of 
cancellation being dependent on the mapping described 
by T. 

Now consider interferers that are spatial point sources 
but emit broadband energy on the band (o^ ^ w ^ oit>* 
The response of the w., branch to an interferer at 0\ Is 
w^6{$^,ut) " gi(bi). To achieve total cancellation Wm must 
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be rhosen lo sairsfy 

w*CsV*ti{H-,,to) - gi(<ij) io^ rs a, ^ (ot,. (6.3a) 

Define the response of each column oi T as 

Uita) - [T]!MfOi. oi) 1 r?: i s M {6.3b) 

where (Tj. denotes the ith column of T, (6.3a) requires 
gi{<u} to be expressed as a linear combination uf the i,(<o], 
1 s£ i :£ M. on ki^ ru £ atty. In general, this cannot be 
accomplished and we conclude that total cancellation of 
broadband interference cannot be obtained. The output 
power due to the broadband interferer can be expressed 
as the integral over frequency of the magnitude squared of 
the diftereiire between the Wq branch and adaptive 
branch resoonscs weighted by the interferer power spec- 
trum. The degree of cancellation can vary dramatically and 
is critically dependent on the interferer direction, fre- 
quency coniont, and choice for T Good cancellation can 
be obtained in some situations when M = 1, while in oth- 
ers even large values of M result in poor cancellation. 
These conclusions are also valid for narrowband sources 
Hidl are broad in direction (spatially distributed radiation). 

B. Partially Adaptive Bcamformer Design. 

The preceding discussion indicates that the degree of 
inierfercncc cancellation is critically dependent on the 
ability of the adaptive channel (Twm) to match the main 
beam response over the interferer frequency extent. This 
provides a means by which to evaluate partially adaptive 
heamrormers. The majority ol work reported on partially 
adaptive beamforming has been <:oncerned with narrow- 
band environmenls. We begin with a discussion of severaf 
narrowband approaches and briefly discuss their ex- 
tension to broadband situations. We then consider 
techniques that are directly applicable to broadband situ- 
ations. Some techniques select for a given C„ while 
others select T directly.^ 

Several approaches lo reducing degrees of freedom are 
based on processing a subset of the outputs of the matrix 

C. ,. Ihis implies that Tv. is a sparse matrix of zeros and 
ones. The outputs ot C. in the MSC are simply auxiliary 
sensor outputs. Morgan 1 1978] evaluated partrally adaptive 
beamformer performance when Tm selected various sub- 
sets of the auxiliary outputs. This is termed an element 
space approach, since a subset of the sensor element out- 
puts is utilized. Several investigators, including Vural 
119771, Adams, el al. [19801. and Gabriel [1986al have con- 
sidered choosing the columns of T to form beams. This is 
traditionally termed a beam space approach. The columns 
of C„ are designed as data independent beamformers, 
each steered to a different direction, and Tm can be used 
to select a subset of the beam outputs. The objective is to 
direct a Ijeam at each Interfering source so that it can be 
subtracted from the output of the w„ branch. One way lo 

■■In priactplr or\o can generate auxiliary constraint!; in on LCMV 
bramrtirmer to reduce the number of adaptive vv»fightb in a CSC 
implementation (Crittiths (19871). Here we assume all mnstraiots are 
alrftady .spcc.iiied in partialJy adaptive LCMV beamlVirming, 



accomplish this is by selecting enough beams to cover all 
possible directions from which interferers might arrive. 
Another is to utilii:e source direction finding techniques to 
select which beams correspond to estimated interferer 
directions. The biggest advantage of the element space 
approach is the simplicity of implementation. Improved 
performance obtained using beam space processing is es- 
pecially evident for interference due to either spatially 
distributed sources or sources with appreciable temporal 
bandwidth. However, this improvement is obtained at the 
expense of impfementing the required beams. 

Chapman I1976J and Owsley 11978] have considered 
choosing the columns of T to select subarrays> i.e... each 
column involves only a subset of the sensors In the array. 
The weightings applied to* each subarray (elements of T) 
can be chosen in various ways, one of which is to use the 
subarray to form a beam. Performance depends on the 
number of sensors in each subarray. which sensors are 
used in each subarray, and the weightings used to com- 
bine the sensor outputs in each subarray. Note that each 
column of T will have zeros in locations corresponding to 
sensors excluded from that subarray, so the overall T is of 
.sparse structure. 

Owsley 11985] suggests a narrowband mcthud for the 
CSC in which the columns of Tm are chosen as a basis for 
the space spanned by the fully adaptive wei]u;ht vectors. 
The dimension of (his space is given by the rank of the 
spatially correlated component of the interference co- 
variance matrix. Van Veen 11988a] extended this approach 
to the broadband case. The dimension of the fully adap- 
tive space can be large in this case since it is given by the 
rank of the correlated components of the broadband in- 
terference covariancc matrix. 

The above approaches are capable of satisfactory per- 
formance wjth narrowband interference since cancella- 
tion requires about one degree of freedom per interferer 
to match the main beam response at each interferer direc- 
tion. However, with broadband interference the main 
beam response must be matched over a range oi fre- 
quency at each interferer direction. Althoujfh the narrow- 
band approaches can be extended, it is difficult to do this 
and keep the number of adaptive weights M small. For 
example, several banks of beams could be designed to 
span the range of directions^ with each bank operating at 
different frequency. The problem is that the number of 
beams required can become large as the frequency band- 
width increases. We seek a T that is efficient, i.e., provides 
good cancellation with a minimum of columns. 

Van Veen and Roberts n9fl7al have considered an opti- 
mization based approach for choosing the columns of T in 
the context of LCMV beamforming. The matrix C„ is de- 
signed to meet the constraints, reducing the problem to 
an unconstrained optimization over the elements of Tm. 
Tm is chosen to minimize the average interference output 
power over a range of likely interference environments. 
Let a vector or parameterize the interference environment 
of interest. In general « can represent interferer loca- 
tions, power levels, spectral distributions, numbers of in- 
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terferers, etc. Defining P](a) as the interference output 
power, Tm »s chosen according to 



P,(a)d£jr 



(6.4) 



where [a^, ut>] denotes the set of interference? scenarios of 
interest. Since output power corresponds to the error be- 
Kvecn the branch and adaptive channe! responses, in 
effect (6.4) selects Tm to provide the best response match 
possible for interference environments in {a»,otuh 

Equation (fa.4) represents a design problem that is 
nonlinear in the design parameters (elements of Tm). A 
suboptimal approach to solving (6.4), in which Tv^ is se- 
quentially designed one column at a time, has been shown 
to be effective in achieving near fully adaptive interference 
cancellation using a small adaptive dimension. The prob- 
lem is still nonlinear; however an effective approximate 
solution is obtained by solving a linear least squares prob- 
lem. An interpretation of this approximate solution is 
given in Van Veen 11988b]. 

Fig. 6.1 illustrates the magnitude response at several 
frequencies for a partialiy adaptive beamformer designed 
using this approach tor the same interference environ- 
ment and array geometry as the fully adaptive example 
in Section IV fsee Fig. 4.2). The partially adaptive beam- 
former has 28 adaptive weights compared to 70 for the 
fully adaptive beamformer. Fig. 6.2 compares the magni- 
tude response of the fully and partially adaptive beam- 
formers as a function of frequency at the interterer DOA's. 
Tm is designed assuming « Is a two-dimensional vector 
with each element corresponding to the DOA of an inler- 
ferer. Both intcrferers have while Specira on 27r/5 ^ 
An/r* And power levels of 40 dB relative to white noise. The 




Figure 6. 1; This.plDt <iepicts the magnitude response of 
a; pdr^iatly adaptm-beamformBr designed to minimize av- 
eragei output power at eight fivquencies on the normal- 
ized fi^uency interval t2«ir/5.4ir/5]. The beamformer 
.uses; ^26: adaptive weights and was designed Bssuming 
tiwvb 'Interfeners. of pow^ 40 dB relative to white noise 
arlvin 9: between -45 and 45 (te^ea The actual inter- 
ferenca/enviranmant. constraint 5, and array geometry 
are ttio same:aa in Fig. 4.2. Note that the beamformer of 
Flg;.4.2 us«s=70 adaptive weights. 
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Figure 6.2 This plot depicts the magnitude response of 
the partially and fuHy adaptivB.beamf ormare correspond- 
ing to Figs. 4.a and 6.1 as a functiian of frequency at £ha 
interferer DDA's. 



region [a,,ab] includes all DOA's between -45 and 
45 degrees. Partially and fully adaptive array gain was 
evaluated at 190 points in this region; the maximum dif- 
ference between partially adaptive and fully adaptive 
array gain was 2.4 dB and the average .7 dB for this par- 
ticular design. 

VII. BEAMFORMER IMPLEMENTATIONS 

In ils simplest form, a beamformer represents a linear 
combination of the sensor data. However, there are many 
different ways of implementing the weighted sum, each 
with its own performance and complexity characteristics. 
The section begins by introducing a general description 
for implementations and then briefly overviews several 
beamformer implementations. For additional detail the 
reader is encouraged to consult the references given in 
this section. 

In essence, most implementations can be represented 
by decomposing w into a product of matrices and a vector 



nv. 



(7.1) 



Vj, 1 ^ i ^ Vo, are a series of matrix transformations of 
conformable dimensions and Wv Is a vector. As a general 
rule, the matrix transformations are chosen to enhance 
performance and/or lower computalioruil complexity. The 
FFT implementation of the DFT Is analogous to (7.1), since 
the DFT matrix can be expressed as a product of permuta- 
tion and butterfly matrices, representing a series of very 
simple computations. The CSC is an example of (7.11 where 
Vo = 1, Vn = Iw«|Cnl, and Wv = ri/Wj^r'. Vi simplifies 
adaptive algorithm implementation by permitting uncon- 
strained adaptation of Wn- In order to simplify multiplica- 
tion of the data by Wy and Cn, several researchers have 
studied decomposition of (WojC J into products of simple 
matrices (e.g., Kafson and Yao 119651 and Ward, et al, [19B6J). 
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As discussed in SccUon II, broadband beam forming can 
be pLTiormcd in the froquenc>' domain or lime domain. 
The vvoighis used lo obtain the outputs at each frequency 
for the frequency domain Ijcamformer depicted In Fig. 2.4 
are easily represented in Terms of (7.1) by employing the 
matrix representation of the DFT. Fri:qucncy domain opti- 
mum bedmformors usually choose rhe weights at each 
frequency based solely on the data at the frequency. This 
partitioning of the data influences both performance and 
computational complexity. Discussion and evaluation of 
frequency domain bejmforming is given in Owsley [198.')], 
Vural .19771, and Gabriel [1986b|, 

Sysiolic implementations of optimum beamformcrs 
have been studied by a number of investigators- They arc 
usually designed to both compute and implement the 
adaptive weights. In general., systolic implementations 
can be des<:ribed in terms of (7.1) where each V has a 
structure amenable lo paratlcl computation and local com- 
munication. MrWhirter (iya3I (see also Haykin [1986] 
<:h. lOi has developed a systolic array that compute? the 
bnamformer output vvithr>ut explicit computation oi the 
adaptive weight vector. Additional references to systolic 
Implementations Include Schreiber and Kuekes [I985I. 
Ward, ct al. !19Hf>l, Owsley 1 1987 1, and Van Veen and 
Roberts |1987b|. 

The study of beamformer implementations is an evolv- 
ing research area, future developments will result from 
advances in VLSI and parallel computing technologies. 

VIII. SUMMARY 

A beamformer forms a scalar output signal as a weighted 
combination of the data received al an array of sensors. 
The weights determine the spatial filtering characierisiics 
of tlie beamformer and enable separation of signals having 
overlapping frequency content if they originate from 
different locations. The weights in a data independent 
beamformer are chosen lo provide a fixed response inde- 
pendent of the received data. Statistically optimum beam* 
formers select the weights to optimize the beamformer 
response based on the statistics of the data. The data sta- 
tistics arc often unknown and may change with time so 
adaptive algorithms are used to obtain weights that con- 
verge to the statistically optimum solution. Computational 
considerations dictate the use of partially adaptive beam- 
formers with arrays composed of large numbers of sen- 
.sors. Many different approaches have been proposed for 
implementing optimum beamformcrs. Future work will 
likely address signal cancellation problems, further reduc- 
tions in computational load for large arrays and improved 
structures for implementation. Beamforming truly repre- 
sents a versatile approach to spatial filtering. 
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